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ABSTRACT 

The mismatch in coefficients of thermal expansion (CTE) of fiber and matrix in metal 
matrix composites reinforced with ceramics fibers induces high thermal stresses in the 
matrix. Elasto-plasdc analyses - with different degrees of simplification and modelization - 
show that an interface layer with a sufficiently high CTE can reduce the tensile hoop stress 
in the matrix substantially. 


INTRODUCTION 

Metal matrix composites reinforced with ceramics fibers have attractive properties for 
engineering applications. These composites have a high strength associated to a low 
density, the ceramics fibers provide a high temperature resistance, and their brittle behavior 
is compensated to some extent by the ductility of the metal matrix. However, because 
ceramics have a low coefficient of thermal expansion (CTE) and metals have a higher CTE, 
the induced thermal mismatch is responsible for residual stresses in the composite when 
subjected to a change in temperature. 

Cracking in the matrix has been observed after cooling down from processing 
temperature to room temperature for brittle matrix materials [1]. It has been proposed that 
the insertion of an adequate interface layer between the fiber and the matrix can reduce the 
tensile stresses in the matrix to a level which avoids matrix cracking. 

Some numerical parametric studies ([2], [3]) suggested that the optimum interface 
layer should have a CTE between those of the matrix and fiber, with a low Young's 
modulus and a high thickness. However, these conclusions were based on a questionable 
optimization procedure [2] or on limited numerical results and did not provide an 
understanding of the problem. 

Recently, Jansson and Leckie [4] conducted a simplified analysis assuming a rigid 
fiber and a very thin layer. They concluded that a compensating layer with a sufficiently 
high CTE can reduce significantly the residual stresses in the matrix. A complete elastic 
analysis was performed by Doghri et al [5], who proposed an optimization procedure 
offering a window of candidate layer materials. Both of the two studies showed that, while 
.the hoop stress in the matrix can be reduced substantially, the axial stress in the matrix is 
less affected by a layer. They also pointed out that plastic yielding may occur in the 
interface layer. 

In this paper, plasticity is taken into account. The first pan of the study is based on a 
three cylinder model, isolating one fiber with an interface layer and a matrix layer (Fig. 1). 
Only thermal loading is considered. The interface layer is elastic-perfecdy plastic while the 


matrix is considered either elastic or elasto-plasdc. A remarkably accurate and efficient 
numerical procedure is developed allowing to have approximations of the stresses and 
strains fields everywhere in the composite. 

In the second pan of this study, the accuracy of the three cylinder model is 
investigated by studying the response of a unit cell of a hexagonal array of the fibrous 
composite. Finite element computations are performed under seveal thermal-mechanical 
loadings. 
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NOTATIONS 


f fiber 

l layer 

m matrix 


j f, l, m 

Rj external radius of j 

t^ = R^-Rf layer thickness 

The radii Rf and R m are related to the fiber volume fraction Cf by Cf 



Ej.Vj 


Young's modulus, Poisson's ratio of j 


Xj, |ij Lam£ coefficients of j 




In 


Hj = 


12 . 


(1 - 2vj)(l + Vj) J 2(1 + vp 
CTE of j 

change in temperature 
radial stress (<%) in j 
hoop stress (Oqq) in j 
axial stress (cr 2Z ) in j 
In the same manner, we define the total strains: Ejj, eqj, e 2 j 

and the plastic strains in the layer 
Oj v. Mises equivalent stress in j 

In the same way that the v. Mises equivalent stress is related to the distorsion energy, a so- 
called damage equivalent stress [6] is related to the total elastic energy and is defined by 


a j 

AT 

a rj 

a 0j 

a zj 


a J=°j 


-a+vp+sa-zvj) 


V 2 

V a j J 


-il/2 


where <?Hj = (<*ij + <*0j + <*zj)/3. 
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1. Development of the Stress-Strain Analysis for the A xisymmetric Thermal Problem 
The composite is not subjected to transverse loading and the outer surface of the 
compound cylinder is traction free. The thermal loading is assumed to be axisymmetric 
with respect to the z-axis, so that the displacement in the transverse plane is radial: Uj(r). 
The fibers are long and the strain and stress distributions are uniform in the z-direction 
except at the end regions which are not studied here. A generalized plane strain assumption 
is made so that e^j * constant = e z . The other strains are given by: 

U; 

e ij = Uj and £qj = — ± 


and there are no shear strains. 

One may use the finite element method (FEM) to solve this problem. But the problem 
begin axisymmetric, we developed a procedure which proved to be more flexible and much 
less computer time consuming than the FEM. The interface layer is assumed to be elasto- 
plastic while the fiber and the matrix are elastic. 

We are looking for a solution which will satisfy the following conditions: 

( 1 ) The constitutive equations everywhere 

(2) Exactly the equilibrium in (0 and (m) 

(3) The continuity of the displacement between (f) and (t ) Uf = at r = Rf 

(4) The continuity of the displacement between (f) and (m) U/ = U m at r = 

(5) The stress free condition at the outer surface of the copound cylinder = 0 at 
r = Rm 

(6) A weak form of the equilibrium in ( t ) 

(7) The axial equilibrium condition associated with the generalized plane strain 



0 
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This solution will not satisfy the continuity of the stress vector at the interfaces and 
(/)/(ra). However, it was found that the resulting jumps are not important. The iterative 
procedure adopted is the following: 

(i) Propose a temperature increment AT 

(ii) For the given temperature increment, propose a total strain increment in 

(t) 

- Find stresses in (/) which satisfy the constitutive equations 

- Find strains and stresses in (0 and (m) which satisfy the conditions 

(1), (2), (3), (4) and (5) 

- Check if the equilibrium equations (6) and (7) are satisfied: 

• If NO, propose another strain increment in (/), i.e. go to (ii) 

• If YES, for this temperature, the solution was found. 

1.1 Equilibrium Conditions (6) and (7) 

Any form of the total strain field in (£) may be assumed, provided that the 
compatibility conditions are satisfied. For example, the following form of the total strains 
in (/) is assumed 


e r/ = c H~^ D r 


f R.Y 


e er = c / + l 


D, 


V r ) 


where c/ and D; are constants to be found. 

We have used the present form of the total strain field in (t) because it has a uniform 
trace, hence the stress field has a uniform trace also and since = constant = e z , the stress 

&zt are likely to be almost uniform, which is in agreement with a generalized plane strain 
assumption. 
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The radial displacement is given as 


U/ = c / r + D / -^ 


The local equations of equilibrium in the layer are 


dr 


+ 




We use the same procedure by which the local equations of equilibrium are extended to the 
global ones. We multiply the local equations by U/ and integrate over the domain of (() to 
give: 


da 


dr 


'7( a rr~°e0) 


r; 

f r+D r^ 


rdr = 0, Vc,, VD, 


This is a weak form of the equilibrium because we assume only a particular form of the 
displacement field. 

The previous equation leads to the following system 


R,r 


R,l 




dr = 0 


2 d a^ R;, . 


dr = 0 


By partial integration this system becomes 
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This system represents the equilibrium condition (6). To this system we add the 
equilibrium condition in the z-direction, i.e. equation (7). Since and are constant, 
we obtain: 



We now define the following matrix notations, which are suitable for numerical analysis. 


V 


l 

-(R,/r ) 2 

0 

c r 

E er 

s 

i 

(R,/r ) 2 

0 

D r 

. - 


0 

0 

1 

_ e z. 


= [B]{e) 



The system of 3 equations of equilibrium can be rewritten as: 
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V J 


where the upper index T denotes the transposition. 

The stresses {o^} are implicit function of {e} and the stresses Ojf, o z f, and <J zm 
are explicit functions of {e}. Then, the last system is nonlinear system of 3 equations and 
3 unknowns which are the components of (e). This method looks like a FEM with only 3 
DOE (degrees of freedom), which makes it a particularly efficient tool. 

We introduce the following matrix notations: 



V 


r 

D ^ 

h. 

r 

0 

V 

0 

= 

0 

0 

0 


ze 


0 

0 

z 

e 

l Zj 





L z J 


= [N]{e) 


ii 

■-Of( R f )- 


*a ( R /)' 

rm £' 

0 

£ 

II 

0 


0 


0 




0 


0 



where { Ff } and {F m } designate the external surface forces applied to the layer due to its 
contact with the fiber and the matrix. 

With these notations, the equilibrium equations can be rewritten again as: 
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{««)} = ? [B] T {a / } rdr - R f (NCR r )] T {F f }- R f [N(R,)) T {F m } ♦ {fJ = (0) 

R f 

It is easily shown that Qther forms of the strains in (/) can be assumed, with the 
corresponding expressions of [B] and [N] being substituted into the last equation. 

1.2 Iterative Approximation of the Total Strain Increment in the Lava 
Keeping the previous notations, we have { Ae/ } = [B] { Ac) . 


For the starting iteration, {Ae} = {0} is proposed. 

If the solution found in (0. (0 and (m) does not satisfy the equilibrium 
equations, the Newton method is used in order to propose another approximation for { Ae) : 


{J(e)} + 


a t _ 

ae 


(C e ) = (0) 


where {C e } means a "correction” to {Ae} and 

R 

% = \ [B] T [H][B]rdr + [h + ] 
ae J R f 

where [H] is the "tangent" modulus in each point of ( i ). If perfect plasticity associated 
with a von Mises yield surface is considered, and if the material properties are assumed not 
to change with temperature, then the "continuum” expression of the tensor H during plastic 
yielding is 




® a 


D 

t 
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where is the Hooke's operator in (/) and is the deviatoric stress tensor in (t). 

The details of the computation of the matrix [h+] are given in Appendix A. This 
matrix takes into account the fact that the external Io4ds applied to (i) are not given 
constant, but depend on (e). 

1.3 Stresses in the Laver 

Once { Ae } is given, one knows {Ae/} everywhere in (/). The stresses (and the 
plastic strains) are then computed by integrating the constitutive equations. The method 
used is an elastic predictor/plastic correction one. A fully implicit scheme is adopted and 
the nonlinear equations are solved by a Newton method. It was shown recently that the 
general procedure presented in [7] becomes particularly simple and efficient: no linear 
system is solved iteratively, the corrections being found explicitly [8]. The details are not 
given here. 

2. Numerical Results of the Axisvmmetric Thermal Problem 

The elastic studies performed by Jansson and Leckie [4] and Doghri et al. [5] proved 
that a compensating interface layer with a sufficiently high CTE can reduce significantly the 
residual stresses in the matrix. They also show that the hoop stress in the matrix can be 
reduced substantially but the axial stress in the matrix is less affected by a layer. These 
studies pointed out that the stresses in the layer can be high enough to induce plastic 
yielding. 

In this section, we will study the effect of a plastic layer on the previous elastic 
results. The numerical results given in sections 2.1, 2.2 and 2.3 were found using the 
method developed in Section 1. 

The fiber and the matrix are defined by the following data: 

Fiber: SiC (SCS6) 

Ef = 360 GPa, v f = 0. 17, ctf = 4.9 x 10*6/*C, Rf = 70 ^m 
Matrix: Ti3Al 

Em = 75.2 GPa, v m = 0.25, a m = 1 1.7 x 10- 6 /‘C, a ym = 380 MPa, 
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R m = 1 10|im 


The radii Rf and R m are related to the fiber volume fraction by 




= 40.5% 


A practical optimization procedure is proposed in [5] offering a window of candidate layer 
materials. According to this procedure, Silver, for example, seems to be a good candidate 
since matrix stresses are considerably reduced. 


The propoerties of an Ag layer and the chosen thickness are 

E t = 71 GPa, v* = 0.30, a ( = 25.9 x lO- 6 /*C, c y / = 100 MPa, t( = lO^m 


The temperature dependence of the materials properties will be taken into account in Section 
2.4. 


2.1 Monotonous Cooling 

A monotonous cooling corresponding to a change in temperature of AT = -800‘C was 
performed with the data above. The final values of stresses and strains at the inner radii are 
given in Tables 1 and 2, respectively. 

As in [5] a reference elastic case, corresponding to a 2 cylinder model (fiber and 
matrix, without layer) is used to measure the changes induced by the interface layer. 

Analyzing the results, the following observations can be made. 

- Plasticity in the layer results in a small increase (less than 10%) of the max. Mises 
stress or damage stress in the matrix by comparison to the elastic layer case. This explains 
why the reductions in these stresses are still important: 
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a 

— E- 2 0.50 if (t) is elastic and 
°ref 


-^12- = 0.55 if (/) is plastic 


The reduction in the max. hoop stress in the matrix is even more important: 


8 * 2 . = 0.37 if (0 is elastic and — SSL = 0.22 if (/) is plastic 
°9ref a 0ref 

- When the layer is added, all the stresses in the matrix are reduced, but one notices 

that: 

Without the layer, the max. tensile stress in the matrix is the hoop stress 

( c ^ 

-52L = 0.84 

<°0m J 

With the layer, the axial stress in the matrix exceeds the max hoop stress and 
this trend is accentuated when plasticity in the layer is taken into account: 

^HL = 1.79 if (() is elastic and = 3.46 if ( i ) is plastic 

a 0m °0m 

Since the stresses in (0 are reduced when plastic yielding is considered, there 
must be an increase in stresses elsewhere. We see that the stresses in (0 decrease and that 
<Jrm and CQ m decrease, then it is to be expected that (%„ increases, by comparison with the 
elastic layer situation. 

One observes that vhen the layer is added, the matrix docs not yield. So the 
assumption which was made at the beginning of this part of the study (elastic matrix) is 
legitimate. 

We have mentioned at the beginning of Section 1 that our solution does not 
verify the continuity of the stress vector at the interfaces (f)/(f) and (()/( m). But the 
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following results show that the jumps in o n are not too important. 


fo rf (R f ) = -43. 57 MPa f 0 r ,(R,) = -36.03MPa 
l a rf(Rf) = _36 - 27MPa [ a nn( R P = -30. 41 MPa 

Table 2 gives the strains in (0, (t) and (m). It can be noticed that even when 
(/) is considered elastic only, the strains in (t) are important (due to high thermal strains). 
When we consider plasticity in (/) the trend becomes accentuated (because the stresses in 
(i) are reduced). 

2.2 Sensitivity Study 

Under a monotonous cooling of AT = -800*C, the influence of two layer parameters 
was studied. 

a) The effect of the variation of the layer thickness was considered, all the other 
parameters being kept constant The two cases: (£) elastic and (/) plastic were studied. 

t. 

Figure 2 shows that if — £ 0. 1 1 , the two cases give the same values of the max. 

R f 

t, 

Mises stress in the matrix. If — > 0.1 1 then for a given thickness, the max. Mises stress 

R f 

in the matrix is higher when the layer is plastic. 

Best reductions in the max. Mises stress are: 


ija. s 


'ref 


0.5 when (O is plastic, and it is obtained for 
0.4 when (l) is elastic, and it is obtained for 


— = 0.16 

R f 

t, 

= 0.23 

Rf 


The same conclusions apply for the damage equivalent stress. 
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b) The effect of the variation of the layer yield stress was considered, all the other 
parameters being kept constant. Since Mises stress = yield stress in (/) (perfect plasticity), 
one knows that if G y / increases then the stresses in (f) increase and one expects that this 
may reduce the stresses in the matrix. But Fig. 3 shows that the value of G yl does not 
affect the Mises stress in the matrix. However, this result leads to another interesting 
prediction: if the temperature dependence of properties is taken into account, plasticity in 
the layer will indeed occur sooner (since a y / will be smaller) but the matrix will not yield 
and the results will not change too much. This will be checked in Section 2.4. 


2.3 Thermal Cycling 

To gain a preliminary idea on how the composite behaves under thermal cycling, a 
linear temperature history (ranging between T! = 825*C and T 2 = 25*C) has been 
considered. The results of only 3 cycles will be presented because it was found that the 
composite has almost shaken down after 3 cycles for the given temperature amplitude. A 
more complete study of the shakedown of the composite will be presented in Section 3.3. 
The materials data are those given at the beginning of Section 2. 

The mechanical strain is defined by: 


-MEC 

ij 



OATS 


= 1 if i = j and 
, S ij = 0 if ‘ * j 


Figure 4a shows that during all the loading history the axial stress in the matrix exceeds the 
max hoop stress and is equal to the max Mises stress. Both the axial stress and max hoop 
stress remain tensile during the cycles. It can be seen from Fig. 4b that the axial 
mechanical strain in the matrix also exceeds the max hoop mechanical strain, but their 
values remain small. 

Figure 5a shows the stress-strain response in the tangential direction at the inner 
radius of the interface layer. It was found that the values of the hoop stress and axial stress 
in the layer are almost identical. Figure 5b shows that the axial mechanical strain and the 
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hoop mechanical strain have almost the same values also. 


2.4 Temperature i mpendence of the Materials Propert ies 

In this part of. the study, the temperature dependence of the materials properties is 
taken into account 

The following materials were considered: 

Fiber SCS6, Matrix: T^Al and Layer Ag 

The materials parameters were compiled from Ref. [9] to [131. The fiber is elastic and the 
layer and the matrix are assumed to be elasdc-perfectly plastic. The data are given in Figs. 
6a-d. 

The finite element package ABAQUS [14] was used. The F.E. mesh used for the 3 
(or 2) cyclinder models is given in Fig. 7a. A generalized plane strain assumption is made. 
More details about the boundary conditions and the modelizadon are given in Appendix B. 

A monotonous uniform cooling AT = -800*C was applied in each of the two cases: 2 
cylinder model (reference case) and 3 cylinder model. Let us recall that in the previous 
sections, the reference case was considered to be elastic. The final values of stresses and 
strains at the inner radii of each material are given in Tables 3 and 4. These results show 
that: 

in the reference case (fiber and matrix) the matrix yields. The addition of the 
layer avoids yielding in the matrix. 

the layer does reduce the hoop stress in the matrix, and this reduction is 

important — SUL = 0.23. 
a 9ref 

It is remarkable to notice that the reduction in the hoop stress (with the same value) as 
well as the fact that the axial stress in the matrix exceeds the max hoop stress were 
observed in the previous pan of this study when the temperature dependence of the 
materials properties has not been taken into account One also notices that the values of the 
strains are almost identical to those obtained in the previous pan. We may conclude that, 
from a qualitative point of view, taking into account the plasticity in the layer is more 
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important than the temperature dependence of the materials properties. 

Figure 8a shows that the addition of the interface layer reduces the max Mises in the 
matrix and avoids plasticity to occur in the matrix. Fig: 8b shows that the interface layer 
can reduce dramatically the max hoop stress in the matrix, while the axial stress in the 
matrix is increased: 

= 0.23 and = 1.42 

G 0ref °zref 

3. ThermalrMechanical Computations fo r a Hexagonal Array of the Composite 

The manufacture and testing of composite materials is demanding in time and effort 
so that there is special need for procedures which can help determine the macroscopic 
properties of the materials from the properties of the constituents. Success in this endeavor 
would mean that prediction of composite properties would bypass expensive and time- 
consuming manufacture and test procedures. 

The so-called homogenization procedure is used to predict the macroscopic properties 
of the composite from the properties of the individual constituents. 

The composite under consideration is composed of a hexagonal array of SCS6 fibers 
coated with the interface layer, in a TijAl matrix (Fig. 9). A hexagonal array presents more 
symmetries then a square array. Taking into account the structure periodicity, the problem 
can be reduced, for symmetric loadings, to the study of a single unit cell (ODBA), shown 
in Fig. 9, with the appropriate boundary conditions, as detailed in Appendix B. 

Resuls for the unit cell under thermal loading only will allow to appreciate the validity 
of the 3 cylinder model used earlier. With the unit cell, it will be also possible to study the 
composite behavior under transverse mechanical loading since the 3 cylinder model is 
clearly not adequate for non-axisymmetric problems. The study was conducted using 
ABAQUS [14]. The F.E. mesh for the cell is shown in Fig. 7b. 

The matrix and the interface layer are assumed to be elastic-perfectly plastic, while the 
fiber remains elastic. The materials properties of the 3 materials are temperature dependent. 
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3.1 Monotonous Cooling 

A monotonous uniform cooling (AT = -800*0 was applied to the unit cell. The 
resulting values of stresses and strains were found to be almost identical to those obtained 
f OT 3 cylinder model, and presented in section 2.4. So, for axisymmetric problems, the 
3 cylinder model seems to be a simple but accurate representation. 

3.2 First yield poi nts in the matrix 

Several proportional monotonously increasing thermal-mechanical loadings were 
applied to the unit cell. The mechanical loading is an average stress i corresponding to a 
macroscopic uniaxial stress state in the 1 -direction (Ej i * 0, Ey = 0 if i * 1 or j * 1). The 
initial temperature is T 0 = 25*C, and two initial stress states were considered: 
no initial residual stresses, 
initial residual stresses due to a cooling down from 
T i = 825*CtoT 0 = 25*C 

The points of first yield in the matrix for each case are reported in Fig. 10. The results 
show that the residual stresses offer a reserve of elasticity to the matrix and the composite. 


3.3 Shake down 

Thermal cycling is applied to the composite under a constant macroscopic stress Ej j 
in the 1 -direction. Shake down is obtained when the macroscopic strain keeps the same 
value at the end of each cycle. Limit shake down points are reported in Fig. 10 . 

In excess of the shake down limit values, the composite rachets on each cycle of 
thermal loading. Inside the shake down region, the composite behaves elastically. The 
results of Fig. 10 show that the particular interface layer (Ag) considered as an example in 
this work does not extend the temperature range for shake down in the composite. The 
reason is that the yield stress of this layer is too low. As pointed out in [4], to improve the 
low cycle fatigue of the composite, it is desirable to have a layer with a high yield stress. 
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CONCLUSIONS 

The optimization procedure proposed in [5] on the basis of an elastic study allows to 
have good candidate layer materials. The elasto-plasdc analysis conducted here confirms 
that such interface layers can reduce the hoop stress and the v. Mises equivalent stress in 
the matrix significantly. However, the axial stress in the matrix is less affected by a layer 
(and it may even increase). This implies that these layers could be successful in composites 
where predominatly radial cracking is observed in the matrix. The interface layer should be 
ductile enough to sustain relatively high strains. It should also have a yield stress high 
enough to improve the range of shake down of the composite. 

It appears that taking into account plasticity is much more important than the 
temperature dependence of the materials properties. Comparison between the results given 
by the 3 cylinder model and a unit cell of a hexagonal array shows that for axisymmetric 
problems the concentric cylinder is a simple but accurate model. 
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APPENDIX A . Stresses in (f) and (m) and contributions of (f) and (m) to the incremental 
stiffness matrix. 

The fiber and the matrix being assumed here to be elastic, their strain and stress fields 
are completely defined by the knowledge of 3 constants: cf, c m , D n (see [5]). 

The solution in (f) is found by solving eqn. (3): 


c f Rf-CjRf + D, 


Cf = C,+ 


D, 




The radial and longitudinal stresses in (f) at r = Rf are then: 


a rf = 2( *f + ^f )c f + ^f e z “ (3 ^f + 2li f )a f AT 


°zf = 2 ^f c f + + 2 ^f^ c z ~ ( 3 ^f + 2 l A f) a f AT 


The solution in (m) is found by solving eqn. (4) and (5): 


c R. + — ! a D m = c,R, + D,R, 

mzp m it t t 
( 


C m + 
m 


^ R V 

\ R tJ 


D m “ C / + D * 


That gives 


2(X_ + P m ) c m + - 2u D - (3X_ + 2u )cx AT = 0 

m m m m z m m m ^m m 


X e + 2(X +u )(c, + D,)-(3X + 2p. )a AT 

D = m 2 m i L 1 m 

m 


2\x + 2(\+u ) 
m m 


'rJ 2 

m 

v R *; 
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and 


c = c,+ D, - 
mil 




\ R ej 


m 


Then, the radial and longitudinal stresses in (m) at r = R/ are 


cr = 2(X_ + u )c m + Xe - 2u 
rm m m m z ^m 




D m -(3*+2u )a AT 
m m *^m m 


<J m = 2X c +(X ry , + 2|i m )e, -(3X m + 2(i m )a m AT 
zm m m m ^nrr z v m ^ m 7 m 


Computation of [h + ] 

The matrix [h + ] contains the partial derivatives of the stresses Off, Om^R/), Ojf and 
<*zm respect to the components (Cf, D/ and e z ) of {e} (see Sections 1.1 and 1.2). 

• Fiber contribution 
We have 


3c 

& 


Ul, 


3c 

3D, 


f _ 


( D \ 


v R t; 


and 




= 2(X f 


3c f 3c 




+ M 


Jiff 

3D. 



3c f 


WV C 

= 2(X f + ^ f ) — + \ f 


3c, 3c, 3c, 3c 

— ^ = 2^,^ - 2X, 

f 3c ^ 3D ; ‘ 


3c, 


3D 


l 
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MaQ7 x contribution: 


We have 


do.e 3c, 

— + (X.j. + 2^.^) 


3e 


3D 

5c 


2(X +n ) 

m LJ3L_ 


3D 


2 ^m + 2 (\n + ^n.>l 


( n ^ 
"m 


v R *; 


3D« 


3D 

9c 


J3L 


2 ^m + 2 <\n + M 


m 

\ r d 


3c 
3c ^ 


( R V 3D 3c„ 3c 
SL 

\. R t ) 


m s si — m - -| 

dc t 3D/3e z 


EH 

V R tJ 


3D 

3e 


3a 3c 

— no> = 2(X +M- m )tr m " -2 ^ 

3c - m m 3c, 


m 


( R m V3D 

IP 

v R /; 


3c. 


3a 

"3D 


3a 3c 


^R Y 3D 

rn 


m 


v R <; 


■+x 

3c, m 


3a 3c „ 3a _ 

zm = 2\ m = 

3c, m 3c, 3D, 


= 2\ — + (>- m + 2a ) 

3e m 3e m 

z z 
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APPENDIX B . Modelization and Finite Element Discretization 


(a) Two and Thre e Cylinder Models 

Because of the symmetry of the problem, only a quarter of the cross section is 
considered. The dimensions and the mesh used are shown in Fig. 7a. ABAQUS 
generalized plane strain elements CGPE6 and CGPE5 were used, two extra nodes (EX1 
and EX2) are needed to impose the generalized plane strain condition. The node EX1 has 
one degree of freedom (DOF) representing the longitudinal displacement which is the 
increase in thickness of the model. The node EX2 takes care of the rotations of the end 
planes, and these are prescribed to be zero. 

The normal displacements of the nodes of the 2 straight edges of the section in Fig. 
7a are imposed to be zero. The two cylinder model (fiber and matrix) is simply obtained by 
saying that the layer material is identical to the matrix one. 

(b) Unit C dl 

The composite is composed of a hexagonal array of SCS6 fibers coated with the 
interface layer, in a T^Al matrix. The transverse section is given in Fig. 9. The following 
approach considered in [15] was kept in this study. 

Edge effects will not be considered in this study, and hence the analysis can be done 
using the theory of homogenization. We consider symmetric loading such that the edges of 
the unit cell as well as the axes 01 and 02, will remain straight during deformation. Hence 
only a quarter of the unit cell need to be considered for analysis, and this is taken to be the 
positive quadrant of the 01,2 axes. 

If this quarter cell (OFEA) is isolated and allowed to deform symmetrically.it will 
undergo displacements as shown in Fig. 1 la. The centerline DCB, about which the region 
is skew-symmetric, will displace as shown. If (Cul, Cu2) are the displacements of C, 
(Pul, Pu2) are the displacements of a point P on CB, and (Qul, Qu2) are the 
displacements of a point Q on DC such that QC = CP, then: 
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Cul = Cu2 = 0, 

Pul +Qul «0, and 
Pu2 + Qu2 = 0 

If this skew- symmetry is taken into account, only the half (ODBA) of the Quarter cell 
need to be analyzed. The deformed shape of (ODBA) is shown in Fig. 1 lb, after giving it 
a rigid body translation so that points in the side OD do not undergo u2 displacements. 

In analyzing the region (ODBA), the boundary conditions to be used are: 

i. Along OD, u2 = 0 

ii. Along OA, ul=Aul 

iii. Along BA, u2 = Au2 

iv. Along DB, for points such as P and Q, 

Pul + Qul =0 

Pu2 + Qu2 = Au2 

The finite element mesh used is given in Fig. 7b. The type of elements is the same as 
in section (a). 

With the boundary conditions given earlier, the key displacements of the F.E. model 
are the ul and u2 displacements of the node A and the u3 displacement of the extra node 
EX1. In terms of these nodal displacements, the average (macroscopic) direct strains in the 
composite are given by: 

_ ul(A) _ u2(A) _ u3(EXl) 

£ U-~V3^' E 22-— jP' E »~ 1 

2 " 

a unit longitudinal thickness being considered. 
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Stresses in Matrix 


Stress 


IT 

60 

zz 

Mises 

o 

Damage 

o* 

Ref. 

a (MPa) 

-187.3 

442.2 

372.0 

597.5 

602.5 

(0 

o (MPa) 

-50.4 

163.5 

292.8 

300.2 

320.3 

elastic 

o/o ref 

0.27 

0.37 

0.79 

0.50 

0.53 

(0 

a (MPa) 

-30.4 

98.7 

342.0 

327.5 

342.7 

plastic 

o/Crtf 

0.16 

0.22 

0.92 

0.55 

0.57 


Stresses in Layer 


Stress 


IT 

00 

zz 

Mises 

Damage 

(/) elastic 

a (MPa) 

-260.1 

1529.5 

1437.1 

1745.2 

1901.7 

(/) plastic 

c (MPa) 

-36.3 

64.9 

62.5 

100.0 

98.8 


Stresses in Fiber 


Stress 


rr 

00 

zz 

Mises 

Damage 

Ref. 

a (MPa) 

-187.3 

-187.3 

-546.5 

359.3 

536.0 

(£) elastic 

a (MPa) 

-260.1 

-260.1 

-780.5 

520.4 

763.8 

{() plasdc 

a (MPa) 

-43.6 

-43.6 

-417.8 

374.2 

406.6 


Table 1. Stresses at the Inner Radii of the Compound Cylinder after Monotonous Cooling 
(Materials properties given at room temperature) 
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Strains in MATRIX 


Strains 


rr 

80 

zz 

Reference 

e(%) 

-1.45 

-0.41 

-0.52 

(0 elastic 

e(%) 

-1.15 

-0.80 

-0.58 

(l) plastic 

e(%) 

-1.12 

-0.90 

-0.50 

Thermal strains: otmAT 

= -0.94% 




Strains in LAYER 





Strains 


rr 

00 

zz 

(f ) elastic 

e(%) 

-3.69 

-0.41 

-0.58 


e(%) 

-5.27 

-0.38 

-0.50 

(t) plastic 

eP (%) 

3.10 

1.60 

1.50 

Thermal strains: a^AT : 

= -2.07% 




Strains in FIBER 





Strains 


rr 

00 

zz 

Reference 

e(%) 

-0.41 

-0.41 

-0.53 

(/) elastic 

e(%) 

-0.41 

-0.41 

-0.58 

(l) plastic 

e(%) 

-0.38 

i 

O 

OJ 

oo 

-0.50 


Thermal strains: ctfAT = -0.39% 


Table 2. Strains at the Inner Radii of the Compound Cylinder after Monotonous Cooling 
(Materials properties given at room temperature) 
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Stresses in Matrix 
Stress 

rr 

00 

zz 

Mises 

Damage 

Reference a (MPa) 

-130.8 

261.7 

235.4 

a 

380.0 

o* 

377.7 

3 cylinder <J (MPa) 

-13.43 

61.82 

335.0 

317.5 

329.3 

model o/Oref 

0.10 

0.23 

1.42 

0.83 

0.87 

Stresses in Layer 

Stress 

rr 

00 

zz 

Mises 

Damage 

3 cylinder c (MPa) 

-1.40 

98.13 

99.07 

100.0 

117.4 

Stresses in Fiber 

Stress 

rr 

00 

zz 

Mises 

Damage 

Reference o (MPa) 

-159.7 

-159.7 

-503.7 

344.0 

491.3 

3 cylinder o (MPa) 

-41.18 

-41.18 

-426.7 

385.5 

415.8 

model a/Orcf 

0.26 

0.26 

0.85 

1.12 

0.84 


Table 3. Stresses at the Inner Radii of the Compound Cylinder after Monotonous Cooling 
(temperature dependence of materials properties taken into account) 
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Strains in MATRIX 


Strain (%) 

e rr 

CD 
C U 

EZZ 

Accumulated 
plastic strain 

Reference 

-1.61 

-0.437 

-0.517 

0.34 

3 cylinder 

-1.08 

-0.961 

-0.506 

0 

Strains in LAYER 
Strain (%) 

e rr 

e 90 

e zz 

Accumulated 
plastic strain 

3 cylinder 

-5.00 

-0.594 

-0.506 

2.94 

Strains in FIBER 

Strain (%) 

% 

£ 00 

e zz 


Reference 

-0.405 

-0.405 

-0.517 


3 cylinder 

-0.381 

-0.381 

-0.506 



Table 4. Strains at the Inner Radii of the Compound Cylinder after Monotonous Cooling 
(temperature dependence of materials properties taken into account) 


Mises stress/Ref.Mises stress 


3 


2 


0 

0 .05 .1 .15 2 

layer thickness/fiber radius 


.25 


Figure 2. — Influence of the layer thickness on the max Mises stress in the matrix. 
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la ym v ' 


fiber 


matrix 



200 400 600 800 1000 1200 1400 1600 

layer yield stress (MPa) 

Figure 3— Influence of the layer yield stress on the max Mises stress in the matrix. 


STRESSES IN THE MATRIX (MPa) 



(b) Mechanical strains. 

Figure 4.— Thermal cycling - history at the inner radius of the 
matrix cylinder. 












(a) Three cylinder model. 



(b) Unit cell of a hexagonal array. 


Figure 7. — Finite element meshes. 
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104.8 Jim 



£ 1 1 / G ym 


Figure 10.— First yield points in the matrix and limit points of 
shake down of the composite. 


(b) Region of analysis. 
Figure 9.— Fibrous composite. 


Jnit cell parameter: h 


(a) Region of analysis and boundary conditions. 



(b) Unit cell and boundary conditions. 

Figure 1 1— Unit cell definition. 
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